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Abstract 

Causal dynamical triangulations (CDT) can be used as a regularization 
of quantum gravity. In two dimensions the theory can be solved anlytically, 
even before the cut-off is removed and one can study in detail how to take 
the continuum limit. We show how the CDT theory is related to Euclidean 
2d quantum gravity (Liouville quantum gravity), how it can be generalized 
and how this generalized CDT model has a string field theory representation 
as well as a matrix model representation of a new kind, and finally how 
it examplifies the possibility that time in quantum gravity might be the 
stochastic time related to the branching of space into baby universes. 



1 Lectures presented at the 49th Cracow School of Theoretical Physics, "Non- 
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1 



1 introduction 



Presumably quantum gravity makes sense as an ordinary effective quantum 
field theory at low energies. At high energies it is presently unclear how 
to view space-time. Is space-time an emergent low-energy structure as ad- 
vocated in string theory, does it require a new concept of quantization as 
believed by people working on loop quantum gravity, or is quantum grav- 
ity "just" an almost standard quantum field theory with non-trivial (non- 
perturbative) ultraviolet behavior ? As long as we do not know the answer 
we have an obligation to pursue all avenues. 

The theory of quantum gravity which starts by providing an utraviolet 
regularization in the form of a lattice theory, the lattice link length being 
the (diffeomorphism invariant) UV cut-off, and where in addition the lattice 
respects causality, is denoted CDT (causal dynamical triangulations). It 
is formulated in the spirit of the last of the three approaches mentioned 
above: old "boring" quantum field theory with a non-trivial fixed point [1], 
[2] . It allows a rotation to Euclidean space-time, the action used being 
the Euclidean Regge action for the piecewise linear geometry represented 
by the (now) Euclidean lattice (see [3, 4] for details of the Regge action 
in the CDT approach). The non-perturbative path integral is performed 
by summing over lattices originating from Lorentzian lattices with a causal 
stucture. Like in ordinary lattice field theories we approach the contiuum by 
fine-tuning the bare coupling constants. The rotation to Euclidean space- 
time allows the use of Monte Carlo simulations of the theory and in four- 
dimensional space-time, which for obvious reasons has our main interest, 
there exists a region of coupling constant space where the infrared behavior 
of the universe seen by the computer is that of (Euclidean) de Sitter space- 
time [5, 6] (for a pedagodical review, see [7]). One might think this is a 
trivial result, but quite the contrary: in Euclidean space-time the Einstein 
action is unbounded from below and the de Sitter solution is only a saddle 
point. Thus we are clearly dealing with an "emergent" property of the path 
integral, a genuine non-perturbative effect arising from an interplay between 
the regularization and the path integral measure, and it can only be valid in 
some limited region of the (bare) coupling constant space. This is precisely 
what is observed. One of the main questions to be answered in such a lattice 
theory is what happens when the lattice spacing is taken to zero, i.e. when 
the cut-off is removed. Is it possible at all to remove the cut off, and if so 
what are the ultraviolet properties of the theory ? Non-trivial UV properties 
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have been observed [8], properties which have been reproduced by other 
"field theoretical" approaches to quantum gravity [9, 10]. However, in the 
CDT-lattice approach it has so far been difficult to penetrate into the trans- 
Planckian region. Active research is ongoing to achieve precisely this. 

Numerical simulations are very useful for understanding whether a non- 
perturbatively defined quantum field theory has a chance to make sense. 
Likewise it is useful for checking if certain conjectured non-perturbative fea- 
tures of the theory has a chance of being true, and one can even discover new, 
unexpected phenomena. In this way the numerical simulations work like ex- 
periments, and this is the spirit in which the above mentioned simulations 
have been conducted. However, numerical simulations have their limitations 
in the sense that they will never provide a proof of the existence of a theory 
and it might be difficult in detail to follow the way the continuum limit is 
approached since it requires larger and larger lattices. It is thus of interest 
and importance to be able to study this in detail, even if only in a toy model. 
Two-dimensional quantum gravity is such a toy model which has a surprising 
rich structure. Many of the intriguing questions in quantum gravity and in 
lattice quantum gravity are still present in the two-dimensional theory. What 
is nice about the two-dimensional theory is that it can be solved analytically 
even at the discretized level. 

In the rest of this article we will discuss the two-dimensional CDT theory. 
In Section 2 we define the "bare" model. In Section 3 we generalize the model, 
allowing for local causality violations and we discuss the generalized model's 
relation to a specific matrix model (Sec. 4). What is new is that the matrix 
model directly describes the continuum limit of the gravity theory. In Sec. 5 
we show how to formulate a complete string field theory for the model and 
in Sec. 6 we show the equivalence to the matrix model defined in Sec. 4. 
The matrix model allow us to define the theory non-perturbatively and we 
show how to calculate non-perturbatively (i.e. including the summation over 
all topologies) certain observables. In Sec. 7 we discuss relations to other 
models. In Section 8 we show that the string field theory can be understood 
as a special kind of stochastic quantization of space, a phenomenon first 
noticed in the context of 2d Liouville quantum gravity in [11]. Stochastic 
quantization defines a non-perturbative Hamiltonian. This is discussed in 
Sec. 9. 
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Figure 1: Piecewise linear space-time histories 1+1 dimensional quantum 
gravity 

2 The CDT formalism 

CDT stands in the tradition of [12], which advocated that in a gravitational 
path integral with the correct, Lorentzian signature of spacetime one should 
sum over causal geometries only. More specifically, we adopted this idea 
when it became clear that attempts to formulate a Euclidean nonperturbative 
quantum gravity theory run into trouble in space-time dimension d larger 
than two. Here we will discuss the implementation only when the space-time 
dimension is two. 

Thus we start from Lorentzian simplicial space-times with d = 2 and 
insist that only causally well-behaved geometries appear in the (regularized) 
Lorentzian path integral. A crucial property of our explicit construction is 
that each of the configurations allows for a rotation to Euclidean signature as 
mentioned above. We rotate to a Euclidean regime in order to perform the 
sum over geometries (and rotate back again afterward if needed). We stress 
here that although the sum is performed over geometries with Euclidean 
signature, it is different from what one would obtain in a theory of quantum 
gravity based ab initio on Euclidean space-times. The reason is that not 
all Euclidean geometries with a given topology are included in the "causal" 
sum since in general they have no correspondence to a causal Lorentzian 
geometry. 
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Figure 2: The propagation of a spatial slice from time t to time t + 1. The 
ends of the strip should be joined to form a band with topology S 1 x [0, 1]. 

We refer to [3] for a detailed description of how to construct the class of 
piecewise linear geometries used in the Lorentzian path integral in higher di- 
mensions. The most important assumption is the existence of a global proper- 
time foliation. This is illustrated in Fig. 1 in the case of two dimensions. We 
have a sum over two-geometries, "stretching" between two "one-geometries" 
separated a proper time t and constructed from two-dimensional building 
blocks. In Fig. 2 we have shown how to fill the two-dimensional space-time 
between the space (with topology S* 1 ) at time t n and time t n +i = t n + a 
where a denotes the lattice spacing. While we in the lattice model often 
use units where everything is measured in lattice length (i.e. the lattice links 
have length one), we are of course interested in taking the limit a — > to 
recover continuum physics. 

In the path integral we will be summing over all possible ways to connect 
a given Id "triangulation" at time t n and a given Id triangulation at t n+ i to 
a slab of 2d space-time as shown in Fig. 2, and in addition we will sum over 
all Id "triangulations" of S* 1 at times t n . Thus we are demanding that the 
time-slicing is such that the topology of space does not change when space 
"evolves" from time t n to time t n+1 . 

The Einstein-Hilbert action S EH in two dimensions is almost trivial. Ac- 
cording to the Gauss-Bonnett theorem the curvature term is a topological 
invariant and does not contribute to the equations of motion as long as the 
topology of space-time is unchanged. And even if the topology changes the 
change in the curvature term is just a number, the change in the Euler char- 
acteristic of the 2d surface. We will first ignore this term, since we are first 
not allowing topology change. Thus the (Euclidean) action simply consists 
of the cosmological term: 

S| H = A J d 2 x ^ — > Sf^ e = AN 2 (1) 
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where N 2 denotes the total number of triangles in the two-dimensional trian- 
gulation. We denote the discretized action the Regge action since it is a trivial 
example of the natural action for piecewise linear geometries introduced by 
Regge [13]. The dimensionless lattice cosmological coupling constant A will 
be related to the continuum cosmological coupling constant A by an additive 
renormalization: 

A = A + i A a 2 , (2) 

the factor 1/2 being conventional. The path integral or partition function 
for the CDT version of quantum gravity is now 

Gf(Ji,J 2 ;*) = /%]e-« - 
Gf\ Ll ,L 2 ,T) = 

where the summation is over all causal triangulations T of the kind described 
above with a total of T time steps, an "entrance loop" of length li = L x a and 
an "exit loop" of length l 2 = L 2 a. The factor 1/CV is a symmetry factor, 
given by the order of the automorphism group of the triangulation T. 

Our next task is to evaluate the sum over triangulations in (3), if possible, 
analytically. Surprisingly, it can be done [14]. One can simply count the 
number of triangulations in a slice like the one show in Fig. 2 and from this 
the total number of triangulations in T slices. As usual, when it comes to 
counting, it is often convenient to introduce the generating function for the 
numbers one wants to count and first find this function. In our model the 
generating function has a direct physical interpretation. We define 

Gf\x 1 ,X 2 ;t) = e-^e-^GfiL^L^T). (4) 

Thus G ( £\x 1 ,X 2 ; T) is the generating function of the numbers G^(Li, L 2 ; T) 
if we write Z\ = e~ Xl , Z 2 = e~ x ' 2 . But we can also view X as a (bare) dimen- 
sionless boundary cosmological constant, such that a boundary cosmological 
term X ■ L has been added to the action. In this way G ( f > (X 1 , X 2 ; T) rep- 
resents the sum over triangulations where the lengths of the boundaries are 
allowed to fluctuate, the fluctuations controlled by the value X of the bound- 
ary cosmological constant. In general we expect, just based on standard di- 
mensional analysis, the boundary cosmological constants Xi to be subjected 

6 



(3) 



to an additive renormalization when the continuum limit is approached. Like 
(2) we expect 

X = X c + x a, (5) 

where x then denotes the continuum boundary cosmological constant, and 
one, after renormalization, has the continuum boundary cosmological action 
x ■ I. 

We refer to [14] for the explicit combinatorial arguments which allow us 
to find Gf\Xi,X 2 ; T). Let us just state the following results: one can derive 
an exact iterative equation (using notation Z = e~ x , W = e~ Y , Q = e _A ) 

WiZM-^W^M-l) (6) 

This equation can be iterated and the solution written as 

Gf\Z,W;T)= (7) 

F 2t (l - F 2 ) 2 ZW 



[l- ZF) - F 2t+1 (F— Z) ( 1- ZF){\- WF) - F 2T (F- Z) (F— W) 
where F is 

1 - Jl - 4Q 2 

These equations tell us that Q c = 1/2 and that Z c — 1 and we can now take 
the continuum limit in (7) using t = T ■ a and find 

Gi 0) (s,y;t) [A< 



(VX + x) +e- 2 ^( v / A-x) 



X (V\ + x)(V\ + y)-e-' 2 ^ t (V\-x)(VX-yy ^ 
Further, the continuum version of (4): 

POO 

G ( H\x 1 ,x 2 ;t)= / dhdl 2 e-^ h - x ^Gf\h,l 2 ;t). (10) 
Jo 

allows us to obtain: 

cf\i L j^)- e ~ [cothV ~ xt]V ~ Hh+i2) i L ( v^ y (ii) 

sinhvCU h \sinh \/Aw ' 



where h(x) is a modified Bessel function of the first kind. Quite remark- 
able this expression was first obtained using entirely continuum reasoning by 
Nakayama [15]. 

Finally, from (2) and (5) one can now obtain the continuum limit of the 
iteration equation (6): 



d ^(o) , d 



g .G^(x,y;t) + — [^-\)G{ \x,y;t)\=0, (12) 

This is a standard first order partial differential equation which should be 
solved with the boundary condition 

G^(x,y;t = 0) = ^- (13) 
x + y 

corresponding to 

Gf\l 1 ,l 2 ;t = 0)=5(l 1 -l 2 ). (14) 

The solution is thus 

Gf\ x , mt) = ^l^—L^, ( 15 ) 
x 2 — A x(t] x) + y 

where x(t; x) is the solution to the characteristic equation 

^ = -(x 2 - A), x(t = 0) = x. (16) 
It is readily seen that the solution is indeed given by (9). since we obtain 

x{t ) = Vx ^ X + x) '^ t ^ X - x) (17) 
(V\ + x)+e- 2 ^ t (V\-x) 

If we interpret the propagator G^ l 2 ; t) as the matrix element between 
two boundary states of a Hamiltonian evolution in "time" t, 

G^\k,l 2 -t)=<k\e' H "%> (18) 

we can, after an inverse Laplace transformation, read off the functional form 
of the Hamiltonian operator Hq from (12), 

H {x) = ^ {x 2 - A)) , H (l) = -lj^ + Xl. (19) 



This end our short review of basic 2d-CDT. We have here emphasized that 
all continuum results can be obtained by explicit solving the lattice model 
and taking the continuum limit simply by letting the lattice spacing a — > 0. 
The same will be true for the generalized CDT model described below, but 
to make the presentation more streamlined we will drop the explicit route 
via a lattice and work directly in the continuum. 



3 Generalized CDT 

It is natural the ask what happens if the strict requirement of "classical" 
causality on each geometry appearing in the path integral is relaxed. While 
causality is a resonable requirement as an outcome of a sensible physical 
theory, there is no compleling reason to impose it on each individual geometry 
in the path integral, since these are not physical observables. We used it, 
inspired by [12], as a guiding principle for obtaining a path integral which is 
different from the standard Euclidean path integral, which was seemingly a 
necessity in higher than two space-time dimensions. 

In Fig. 3 we show what happens if we allow causality to be violated locally 
by allowing space to split in two at a certain time t, but we never allow the 
"baby" universe which splits off to come back to the "parent" universe. The 
baby universe thus continues its life and is assumed to vanish, shrink to 
nothing, at some later time. We now integrate over all such configurations 
in the path integral. From the point of view of Euclidean space-time we 
are simply integrating over all space-times with the topology of a cylinder. 
However, returning to the original Minkowskian picture it is clear that at 
the point where space splits in two the light-cone is degenerate and one 
is violating causality in the strict local sense that each space-time point 
should have a future and a past light-cone. Similarly, when the baby universe 
"ends" its time evolution the light-cone structure is degenerate. These points 
thus have a diffeomorphism invariant meaning in space-times with Lorentzian 
structure, and it makes sense to associated a coupling constant g s with the 
process of space branching in two disconnected pieces. 

The equation corresponding to Fig. 3 is [16] 

-G x , gs (x,y,t) = -— [(V - A) +2g s W x , 93 (x))g x , 9 M V, 0] ■ (20) 
W\ j9s (x) is denoted the disk amplitude with a fixed cosmological constant x. 
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Figure 3: In all four graphs, the geodesic distance from the final to the 
initial loop is given by t. Differentiating with respect to t leads to eq. (20). 
Shaded parts of graphs represent the full, ^-dependent propagator and disc 
amplitude, and non-shaded parts the CDT propagator. 



It is related to the disk amplitude with a fixed boundary length by 

POO 

W Kga {x) = / dle-* l W x>9s (l)- (21) 
Jo 

It describes the "propagation" of the a spatial universe until it vanishes in 
the vacuum. If we did not allow any spatial branching we would simply have 

poo 1 

Wf (*) = / dt Gf (x, I = 0; t) = 7 =, (22) 

Jo x + V A 

where G^\x,l;t) denotes the Laplace transform of G^(l',l]t) with respect 
to I'. From the composition rules for G\^ gs {li,h]t) it follows that it has 
(mass) dimension 1. Thus Gx^^x^l^t) is dimensionless and it follows that 
the (mass) dimension of the coupling constant g s must be 3. In a discretized 
theory it will appear as the dimensionless combination g s a 3 , a being the 
lattice spacing, and one can show that the creation of more than one baby 
universe at a given time t is suppressed by powers of a (see [16] for details). 
Thus we only need to consider the process shown in Fig. 3. For a fixed 
cosmological constant A and boundary cosmological constants x, y expres- 
sions like G\ t9a (x, y; t) and W\ t9a (x) will have a power series expansion in the 
dimensionless variable 

9s (23) 



A3/2 

and the radius of convergence is of order one. Thus the coupling constant 
g s indeed acts to tame the creation of baby universes and if g s exceeds this 
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Figure 4: Graphical illustration of eq. (66). Shaded parts represent the full disc 
amplitude, unshaded parts the CDT disc amplitude and the CDT propagator. 

critical value eq. (20) breaks down and is replaced by another equation cor- 
responding to Liouville quantum gravity with central change c = 0. This 
was already observed in [14], and elaborated upon in [17] where it was shown 
that one could indeed obtain the CDT from Liouville gravity by integrating 
out baby universes. In fact, in [17] an amazing non-analytic map between 
the two theories was discovered: 



where A and x denote the cosmological and boundary cosmological constants 
in the purely Euclidean theory (Liouville quantum gravity with c = 0). This 
relation was rediscovered in [18] where the CDT parameter xj y/\ was called 
the uniformization parameter of the algebraic surface which in [18] was as- 
sociated to the semiclassical behavior of Liouville quantum gravity. CDT is 
the quantum theory associated with the modified algebraic surface defined 
by the non-analytic mapping (24). Once the coupling constant g s is intro- 
duced the CDT theory is charged, but in an analytic way, until the radius of 
convergence for k is reached. The detailed discussion of what happens when 
k approaches the radius of convergence can be found in [19]. 

Differentiating the integral equation corresponding to Fig. 3 fwith respect 
to the time t one obtains (20). The disc amplitude W\ y g s (x) is at this stage 
unknown. However, one has graphical representation for the disc amplitude 
shown in Fig. 4. It translates into the equation [16] 




(24) 




(25) 



oo oo 



g a J dt J dhdl 2 (h + l2)G{ \x,h + l 2 ;t)Wx, gs (h)Wx, g M 







The superscript (0) to indicate the CDT amplitudes without baby universe 
branching, calculated above. We assume 



W Xt9s=0 (x)=wi°\x), 



(26) 



and similarly for . The integrations in (25) can be performed, yielding 



x + y/X x 2 - X 
Solving for W\ iga (x) we find 

ix 2_ X )+W x , 9s (x) 



(wl gs (V\)-WlJx)). (27) 



W x , 9s (x) = 
where we have defined 



2g s 



(28) 



W x , 9s (x) = y (x 2 - A) 2 + Ag a (g s Wl gs (V\)+x-Vx). 



(29) 



W\ jga (x) is determined up to the value Wx, gs (VX). We will now show that 
this value is fixed by consistency requirements of the quantum geometry. If 
we insert the solution (28) into eq. (20) we obtain 



d ~ 



d 



dt G x , 9s (x,y;t)= ^ 



Wx, gs (x) G Xt g.(x,y;t) 
In analogy with (12) and (15), this is solved by 

W x , 9 Mt,x)) 1 



(30) 



Gx, 9s (x,y,t) 



W x ,Jx) x(t,x)+y : 



(31) 



where x(t,x) is the solution of the characteristic equation for (30), the gen- 
eralization of eq. (16): 



dx 

= -W x ,g s (x), x(Q,x)=x, 



(32) 
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such that 

dy 



(33) 

ix(t) w x ,g s {y) 

Physically, we require that t can take values from to oo, as opposed to just 
in a finite interval. From expression (33) for t this is only possible if the 
polynomial under the square root in the defining equation (29) has a double 
zero, which fixes the function W\ t9a (x) to 

W x , 9s (x) = (x-a) \/(x + a) 2 - 2g s /a = X(x - u) y / ( y x + u) 2 -2K, (34) 
where 

x = xV\, a = uV~\, w 3 — u + k = 0. (35) 

In order to have a physically acceptable W\ jgs (x), one has to choose the 
solution to the third-order equation which is closest to 1 and the above 
statements about the expansion of W\, gs (x) in a power series in k follows 
from (28), (34) and (35). 



4 The matrix model representation 

The formulas (34) and (28) are standard formulas for the resolvent of a 
Hermitean matrix model, calculated to leading order in N, the size of the 
matrix. In fact the following matrix model 

Z{\ 9s ) = j tye-™™, V{<t>) = J s {^~\^) (36) 
has a resolvent 

Qf Tr (^) ) = W XMx) {x) + 0{l/N% (37) 

where W\ t9s (x) is given by (28), and where the expectation value of a matrix 
expression 0(<f>) is defined as 

(0(0) = [ dtp e-^W 0{<t>). (38) 

Z{\,9s) J 

What is surpising here, compared to "old" matrix model approaches to 2d 
Euclidean quantum gravity, is that the large iV limit reproduces directly the 
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continuum theory. No scaling limit has to be taken. The situation is more 
like in the Kontsevich matrix model, which directly describes continuum 
2d gravity aspects. In fact the qubic potential is "almost" like the qubic 
potential in the Kontsevich matrix model, but the wold-sheet interpretation 
is different. 

Can the above correspondance be made systematic in an large N expan- 
sion and can the matrix model representation help us to a non-perturbative 
definition of generalized 2d CDT gravity ? The answer is yes [19]. 

First we have to formulate the CDT model from first principles such that 
we allow for baby universes to join the "parent" universe again, i.e. we have 
to allow for topology changes of the 2d universe, and next we have to check 
if this generalization is correctly captured by the matrix model (36) [20] 



In quantum field theory particles can be created and annihilated if the process 
does not violate any conservation law of the theory. In string field theories 
one operates in the same way with operators which can create and annihi- 
late strings. From the 2d quantum gravity point of view we thus have a 
third- quantization of gravity: one-dimensional universes can be created and 
destroyed. In [21] such a formalism was developed for non-critical strings 
(or 2d Euclidean quantum gravity). In [20] the formalism was applied to 2d 
cdt gravity leading to a string field theory or third quantization for CDT 
which allows us in principle to calculate any amplitude involving creation 
and annihilation of universes. 

Let us briefly review this formalism. The starting point is the assumption 
of a vacuum from which universes can be created. We denote this state |0) 
and define creation and annihilation operators: 



The factor / multiplying the delta-function is introduced for convenience, see 
[20] for a discussion. 

Associated with the spatial universe we have a Hilbert space on the pos- 
itive half-line, and a corresponding scalar product (making H (l) defined in 
eq. (19) hermitian) : 



5 CDT string field theory 





(39) 




(40) 



14 



The introduction of the operators and ^(l) in (39) can be thought of 
as analogous to the standard second quantization in many-body theory. The 
single particle Hamiltonian Hq defined by (19) becomes in our case the "sin- 
gle universe" Hamiltonian. It has eigenfunctions ip n (l) with corresponding 
eigenvalues e n = 2nV% n = 1,2,...: 

V>n(0 = /e-^pn-iW, H (l)ij; n (l) = e„^(Z), (41) 

where p„_i(/) is a polynomial of order n — 1. Note that the disk amplitude 
W^\l), which is obtained from (22), formally corresponds to n = in (41): 

Wi°\l)=e-^ XI , H (1)W<°\1)=0. (42) 

This last equation can be viewed as a kind of Wheeler-deWitt equation if 
we view the disk function as the Hartle-Hawking wave function. However, 
W^\l) does not belong to the spectrum of H (l) since it is not normalizable 
when one uses the measure (40) 

We now introduce creation and annihilation operators a\ and a n cor- 
responding to these states, acting on the Fock- vacuum |0) and satisfying 
K, 4J = Km- We define 

*(J) = $>„^(/), *t( Z ) = ^4<(Z), (43) 

n n 

and from the orthonormality of the eigenfunctions with respect to the mea- 
sure dl/l we recover (39). The "second-quantized" Hamiltonian is 

/°° HI 
j &(l)HoW(l), ( 44 ) 

and the propagator G\{h, h\ t) is now obtained as 

Gi 0) (/ 1 ,/ 2 ;t) = (0|^(/ 2 )e-^'M/t(/ 1 )|0). (45) 

While this is trivial, the advantage of the formalism is that it automati- 
cally takes care of symmetry factors (like in the many-body applications in 
statistical field theory) both when many spatial universes are at play and 
when they are joining and splitting. We can follow [21] and define the fol- 
lowing Hamiltonian, describing the interaction between spatial universes: 

H = H - 9s j dh j dl 2 ^(h)^(l 2 )^(h + l 2 ) (46) 
-ag s J dh j dk&ih + kMl^ih)- J jp(Z)*(Z), 
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Figure 5: Graphical illustration of the various terms in eq. (46). 

where the different terms of the Hamiltonian are illustrated in Fig. 5. Here 
g s is the coupling constant we have already encountered in Sec. 3 of mass 
dimension 3. The factor a is just inserted to be able the identify the action 
of the two <? s -terms in (46) when expanding in powers of g s . We will think of 
a = 1 unless explicitly stated differently. When a = 1 H is hermitian except 
for the presence of the tadpole term. It tells us that universes can vanish, but 
not be created from nothing. The meaning of the two interaction terms is as 
follows: the first term replaces a universe of length l\ + 12 with two universes 
of length l\ and ^2- This is one of the processes shown in Fig. 5. The second 
term represents the opposite process where two spatial universes merge into 
one, i.e. the time-reversed picture. The coupling constant g s clearly appears 
as a kind of string coupling constant: one factor g s for splitting spatial uni- 
verses, one factor g s for merging spatial universes and thus a factor g 2 s when 
the space-time topology changes, but there is also factors for branching alone. 
This is only compatible with an Euclidean sft-picture if we associate a punc- 
ture (and thus a topology change) with the vanishing of a baby universe. As 
discussed above this is indeed not unnatural from a Lorentzian point of view. 
From this point of view the appearance of a tadpole term is more natural in 
the CDT framework than in the original Euclidean framework in [21]. The 
tadpole term is a formal realization of this puncture "process", where the 
light-cone becomes degenerate. 

In principle we can now calculate the process where we start out with m 
spatial universes at time and end with n universes at time t, represented 
as 

G x , gs (h, .., l m ; li, ~, C t) = (0\Wi) ■ ■ -WJ . . . ^(UIO). (47) 

5.1 Dyson- Schwinger equations 

The disk amplitude is one of a set of functions for which it is possible to derive 
Dyson-Schwinger equations (DSE). The disk amplitude is characterized by 
the fact that at t = we have a spatial universe of some length, and at some 
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point it vanishes in the "vacuum" . Let us consider the more general situation 
where a set of spatial universes of some lengths /j exists at time t — and 
where the universes vanish at later times. Define the generating function: 

Z(J) = lim(0|e-^ e f dljmHl) \0}. (48) 

Notice that if the tadpole term had not been present in H Z(J) would triv- 
ially be equal 1. We have 



(49) 



J=0 



Z(J) is the generating functional for universes that disappear in the vacuum. 
We now have 



= lim 



|. ( | e -^ e-W'^lO) = -{0\e- tA HeS dlJ WV\0) 



(50) 



Commuting the \l/(Z)'s in H past the source term effectively replaces these 
operators by U(l), after which they can be moved to the left of any ^(l) 
and outside (0|. After that the remaining \pt(Z)'s in H can be replaced by 
8/8 J (I) and also moved outside (0|, leaving us with a integro-differential 
operator acting on Z(J): 

= ^ dlJ(l)o(l,J,-^J Z(J) (51) 



where 



°{ l - J -jj) = H °Ujk- m (52) 

Z(J) is a generating functional which also includes totally disconnected 
universes which never "interact" with each other. The generating functional 
for connected universes is obtained in the standard way from field theory by 
taking the logarithm of Z(J). Thus we write: 

F(J) = log Z(J), (53) 
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and we have 

}im(0\e- t H%)-..&(l n )\0) cm = , T f/ {J ,\ (] v , (54) 
*->oo 0J(/i) • • ■dJ{l n ) J=Q 

and we can readily transfer the DSE (51)-(52) into an equation for the con- 
nected functional F(J): 

5 2 F(J) 



V = l°°dl J(l) j#o(0 6 ~§^ - - gj J dl' 



5J(l')5J(l - 1') 



. gJ f ^ - ag ,i r *rjw-M& } ■ (55) 



! 5J{V) 5J{l-V) * b J y J 5J(l + l') 

From eq. (55) one obtains the DSE by differentiating (55) after J (I) a number 
of times and then taking J (I) = 0. 

5.2 Application of the DSE 

Let us introduce the notation 

5 n F(J) 



w(h, ...,/„) = 



SJ(h)---SJ(l n ) 



(56) 



j=o 



as well as the Laplace transform w(xi, . . . , x n ) (like in (10)). Let us differen- 
tiate eq. (55) after J (I) one and two times, then take J (I) = and Laplace 
transform the obtained equations. We obtain the following equations (where 
H (x)f(x) = d x [(x 2 -\)f(x)}): 

= H (x)w(x) — 1 + g s d x (w(x, x) + w(x)w(x) j , (57) 



+2g s (d x [w{x)w{x, y)]+d y [w(y)w{x, y)}) + 2ag s d x d y f w ^ x "> ^ w ^ \ 



= (H (x) + H (y))w(x,y)+g s d x w(x,x,y) + g s d y w(x,y,y) (58) 

'w(x) —w(y) ' 
x-y 

The structure of the DSE for an increasing number of arguments is hopefully 
clear (see [20] for details). 

We can solve the DSE iteratively. For this purpose let us introduce the 
expansion of w(x±, . . . , x n ) in terms of the coupling constants g s and a: 



(xi,...,x n )= ^ ® k ^2 9 ™ w(xi,...,x n ;m, k). (59) 



k=n— 1 m=k— 1 
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The amplitude w(xi, . . . , x n ) starts with the power (ag s ) n ~ 1 since we have to 
perform n mergings during the time evolution in order to create a connected 
geometry if we begin with n separated spatial loops. Thus one can find the 
lowest order contribution to w(xi) from (57), use that to find the lowest order 
contribution to w(xi,x 2 ) from (58), etc. Returning to eq. (57) we can use 
the lowest order expression for w(xi,x 2 ) to find the next order correction to 
w(xi), etc. 

As mentioned above the amplitude w(xi, . . . ,x n ) starts with the power 
(ag s ) n_1 coming from merging the n disconnected spatial universes. The 
rest of the powers of ag s will result in a topology change of the resulting, 
connected worldsheet. From an Euclidean point of view it is thus more 
appropriate to reorganize the series as follows 

oo 

w(xi, ...,x n ) = (ags)"' 1 ^2(ag 2 s ) h w h (x u ...,x n ) (60) 

h=0 

oo 

w h {xi, ...,x n ) = y^gfw(xi, . . . ,x n ;n-l+2h+j,n-l+h) (61) 
j=o 

and aim for a topological expansion in ag 2 , at each order solving for all 
possible baby-universe creations which at some point will vanish into the 
vacuum. Thus Wh(xi, . . . ,x n ) will be a function of g s although we do not 
write it explicitly. The DSE allow us to obtain the topological expansion 
iteratively, much the same way we already did as a power expansion in g s . 



6 The matrix model, once again 

Let us consider our N x N Hermitian matrix with the qubic potential (36) 
and define the observable 

W( X1 , .... = «~ ( (Tr (-1^) . . . (tr ) , (62) 

where the subscript d refers to the fact that the correlator will contain 
disconnected parts. We denote the connected part of the correlator by 
W{x\ . . . ,x n ). It is standard matrix model technology to find the matrix 
model DSEs for W(x 1 . . .,x n ). We refer to [22, 23, 24, 25] for details. One 
obtains precisely the same set of coupled equations as (57)-(58) if we identify: 

a = Jpi ( 63 ) 
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and the discussion surrounding the expansion (60) is nothing but the stan- 
dard discussion of the large N expansion of the multi-loop correlators (see 
for instance [25] or the more recent papers [26, 27, 28]). Thus we conclude 
that there is a perturbative agreement between the matrix model (36) and 
the CDT SFT in the sense that perturbatively: 

W(x 1 ,...,x n )=w(x 1 ,...,x n ). (64) 

In practice the SFT is only defined perturbatively, although in principle we 
have available the string field Hamiltonian. However, we can now use the 
matrix model to extract non-pertubative information. The identification of 
the matrix model and the CDT SFT DSEs were based on (63), but in the SFT 
we are interested in a — 1, i.e. formally in TV = 1, in which case the matrix 
integrals reduce to ordinary integrals. This means that we will consider the 
entire sum over topologies "in one go" : 



Z(g, A) = J dm exp 
while the observables (62) can be written as 



1 (^ 1 3 
— Am m 

9s \ 3 



(65) 



exp 



— y s (Am — | m 3 ) 



W ^--^=Z^X)J dm {xi - m) ... (x ,- m) ■ («> 

These integrals should be understood as formal power series in the dimen- 
sionless variable k defined by eq. (23). Any choice of an integration contour 
which makes the integral well defined and reproduces the formal power series 
is a potential nonperturbative definition of these observables. However, dif- 
ferent contours might produce different nonperturbative contributions (i.e. 
which cannot be expanded in powers of t), and there may even be nonper- 
turbative contributions which are not captured by any choice of integration 
contour. As usual in such situations, additional physics input is needed to 
fix these contributions. 

To illustrate the point, let us start by evaluating the partition function 
given in (65). We have to decide on an integration path in the complex plane 
in order to define the integral. One possibility is to take a path along the 
negative axis and then along either the positive or the negative imaginary 
axis. The corresponding integrals are 

Z( 9s ,\) = VXk 1 / 3 F ± (k- 2 / 3 ), F ± (/ t - 2 / 3 ) = 2vre ± ^ 6 Ai(/ t - 2 / 3 e ±2 ^ 3 ), (67) 
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where Ai denotes the Airy function. Both F± have the same asymptotic 
expansion in k, with positive coefficients. Had we chosen the integration 
path entirely along the imaginary axis we would have obtained (27ri times) 
Ai(/t~ 2 / 3 ), but this has an asymptotic expansion in k with coefficients of 
oscillating sign, which is at odds with its interpretation as a probability 
amplitude. In the notation of [29] we have 

F±(z)=n(Bi{z)±iAi(z)}, (68) 

from which one deduces immediately that the functions F±(k~ 2 / 3 ) are not 
real. However, since Bi(/t~ 2 / 3 ) grows like for small k while Ai(/t~ 2 / 3 ) falls 
off like e~3^, their imaginary parts are exponentially small in 1/k compared 
to the real part, and therefore do not contribute to the asymptotic expansion 
in k. An obvious way to define a partition function which is real and shares 
the same asymptotic expansion is by symmetrization, 

I (j F + + F _) = 7rB i. (69) 

The situation parallels the one encountered in the double scaling limit of the 
"old" matrix model [32], and discussed in detail in [34], but is less compli- 
cated. We will return to a discussion of this in the next section. 

Presently, let us collectively denote by F(z) any of the functions F±(z) 
or irBi(z), leading to the tentative identification 

Z{g„\) ^^F^ 3 ), F"{z) = zF(z), (70) 

where we have included the differential equation satisfied by the Airy func- 
tions for later reference. In preparation for the computation of the observ- 
ables Wd(xi, . . . ,x n ) we introduce the dimensionless variables 

x = xV\, m = gl /3 p, W d (x 1 ,...,x n ) = A _n/2 w d (xi, ...,x n ). (71) 

Assuming Xk > 0, we can write 
1 f°° 

We can use this identity to re-express the pole terms in eq. (66) to obtain 
the integral representation 



/•oo n F\K 3— K,lY^ =1 ai) 

W d (%l, ■■■,Xn)= H ^ e-^ iai+ - +S - a ») 7 L 

i\ f(«-§) 



(73) 
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for the amplitude with dimensionless arguments. By an inverse Laplace 
transformation we thus obtain: 



W d {k,...,l n ) - F(k- 2 /3) • ^ 

For the special case n = 1 we find 

Wil) = Ff^PT) (75) 

for the disc amplitude, together with the remarkable relation 

W d (h,...,l n )=W(k + --- + Q. (76) 
By Laplace transformation this formula implies the relation 

W d (x!, ...,*„)=£ U n W{Xl) _ y (77) 
i=1 [[j^i\ X j X 'i) 

From Wd(xi, . . . ,x n ) we can construct the connected multiloop functions 
W(x\, . . . ,x n ) using standard field theory Let us remark that the asymp- 
totic expansion in k of W(x±, . . . , x n ) of course agrees with that obtained by 
recursively solving the CDT Dyson- Schwinger equations. 

7 Relation with other models 

Can we identify a know continuum conformal field theory coupled to 2d 
gravity which leads to the functions we have here calculated? The answer is 
yes. As a starting point note that for the disk amplitude we have 

from which we conclude that the susceptibility exponent 7 = 1/2 (the lowest 
non- analytic power of A in the disk amplitude is A 1-7 ). 7 = 1/2 is the generic 
value of the susceptibility exponent for the so-called branched polymers and 
the way branched polymers enters into the game of 2d gravity coupled to 
conformal field theories is as follows: recall that the (2, 2m — 1), m — 2, 3, . . ., 
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minimal conformal field theories coupled to 2d Euclidean quantum gravity 
can be described as double-scaling limits of one-matrix models with certain 
fine-tuned matrix potentials of order at least m + 1. Formally, the case 
m — 1, which corresponds to a somewhat degenerate (2, 1) conformal field 
theory with central charge c = — 2 (which when coupled to 2d gravity is called 
topological gravity), is then described by a special double-scaling limit of the 
purely Gaussian matrix model (see the review [35]). In this double-scaling 
limit one obtains for the so-called FZZT brane precisely the Airy function, 
see [36, 37, 34, 38] for recent discussions. For this model 7 = — 1. While it is 
possible to describe 2d topological quantum gravity by a double-scaling limit 
of the Gaussian matrix model, the most natural geometric interpretation of 
the Gaussian matrix model is in terms of branched polymers, in the sense 
that the integral 

f dM tr M 2n e~\ trM2 t s 

± (79) 

JdMe~^ M2 V ' 

can be thought of as the gluing of a boundary of length n into a double- 
line branched polymer of length n. Since the branched polymers are also 
allowed to form closed loops, their partition function contains a sum over 
topologies 'en miniature', and one can indeed define a double-scaling limit 
of the model. When solving for the partition function in this limit, one 
obtains precisely our Z(g s ,X) of eq. (70)! (see [39], where this remarkable 
result was first proved, for details). It does not imply that the generalized 
CDT model is just branched polymers, it is much richer since it has a 2d 
surface representation and many more observables, but the non-perturbative 
branching process is clearly that of branched polymers. The relation between 
the generalized CDT model and the c = —2 model is as follows: we have the 
"right" branch of Liouville theory with susceptibility exponent 7_ (= -1) and 
a "wrong" branch with susceptibility exponent 7+, related by 

7+ = - T ^-. (80) 

The interpretation of this 7+ in relation to 7_ in terms geometry can be 
found in [40, 41, 42], and for earlier related work see [43, 44]. It is all 
related to dominance or non-dominance of branching of baby universes. The 
simplest example is precisely given by c = —2: topological quantum gravity 
has 7_ = — 1 whose dual is the "wrong" 7+ = 1/2, which happens to be 
the value occurring generically in the theory of branched polymers (see, e.g., 
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[45, 46, 47] for a discussion of why branched polymers and baby universes are 
generic and even dominant in many situations in non-critical string theory 
and even in higher dimensional quantum gravity). 



8 Stochastic quantization 

It is a most remarkable fact that the above mentioned result can all be 
understood as a result of stochastic quantization of space. In this picture 
time becomes the stochastic time related with the branching of space into 
baby universes and the original CDT model described in Sec. 2 becomes the 
classical limit where no stochastic processes are present [48]. 

Recall the Langevin stochastic differential equation for a single variable 
x (see, for example, [49, 50]). 

i»(t) = -f(x(»)( t )}+Vnis(t), (si) 

where the dot denotes differentiation with respect to stochastic time t, u{t) is 
a Gaussian white- noise term of unit width and f(x) a dissipative drift force: 

m - ^ (82) 

The noise term creates a probability distribution of x(t), reflecting the as- 
sumed stochastic nature of the noise term, with an associated probability 
distribution 

P(x,x ;t) = (6(x-xM(t;x ))) v , (83) 

where the expectation value refers to an average over the Gaussian noise. 
P(x, x ; t) satisfies the Fokker-Planck equation 

= £ + /<*)'<*. *«o). W 

This is an imaginary-time Schrodinger equation, with y/Ti playing a role 
similar to h. It enables us to write P as a propagator for a Hamiltonian 
operator H, 

P(x,x ;t) = (x\e- t6 \x Q ), H =^lf + ipf(x), (85) 
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with initial condition x(t = 0) = x , and p = —id x . It follows that by 
defining 

G(x Q ,x;t) = ^-P(x,x ;t) (86) 

OXq 

the function G(x , x; t) satisfies the differential equation 
dG(x ,x;t) d I 1 dG(x , x; t) 



dt dx \ 2 dx 



f(x )G(x ,x;t) . (87) 



Omitting the noise term corresponds to taking the limit Q — > 0. One can 
then drop the functional average over the noise in (83) to obtain 

d 

P d (x,x ;t) = S(x - x(t,x )), G d (x ,x;t) = — S(x - x(t,x )). (88) 

OXq 

It is readily verified that these functions satisfy eqs. (84) and (87) with Q = 0. 
Thus we have for S(x) = —\x + x 3 /3: 



dG d (x ,x;t) d 



dt dx 



((\-x 2 )G cl (xo,x;t)y (89) 



Comparing now eqs. (12) and (89), we see that we can formally re- 
interpret G^\xo,x;t) - an amplitude obtained by nonperturbatively quan- 
tizing Lorentzian pure gravity in two dimensions - as the " classical prob- 
ability" G c i(xo,x;t) corresponding to the action S(x) = —\x + x 3 /3 of a 
zero-dimensional system in the context of stochastic quantization, only is 
the boundary condition different, since in the case of CDT x is not an or- 
dinary real variable, but the cosmological constant. The correct boundary 
conditions are thus the ones stated in eqs. (13)-(14) 

Stochastic quantization of the system amounts to replacing 

G { °\x ,x;t) ^G{x ,x;t), (90) 

where G(x ,x;t) satisfies the differential equation corresponding to eq. (87), 
namely, 

We have introduced the parameter g s := 0/2, which will allow us to repro- 
duce the matrix model and SFT results reported above. 
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A neat geometric interpretation of how stochastic quantization can cap- 
ture topologically nontrivial amplitudes has been given in [51]. Applied to 
the present case, we can view the propagation in stochastic time t for a given 
noise term u{t) as classical in the sense that solving the Langevin equation 
(81) for x^\t) iteratively gives precisely the tree diagrams with one exter- 
nal leg corresponding to the action S(x) (and including the derivative term 
x^(t)), with the noise term acting as a source term. Performing the func- 
tional integration over the Gaussian noise term corresponds to integrating 
out the sources and creating loops, or, if we have several independent trees, 
to merging these trees and creating diagrams with several external legs. If 
the dynamics of the quantum states of the spatial universe takes place via the 
strictly causal CDT-propagator G = e~ tH ° , a single spatial universe of length 
/ cannot split into two spatial universes. Similarly, no two spatial universes 
are allowed to merge as a function of stochastic time. However, introduc- 
ing the noise term and subsequently performing a functional integration over 
it makes these processes possible. This explains how the stochastic quanti- 
zation can automatically generate the amplitudes which are introduced by 
hand in a string field theory, be it of Euclidean character as described in [51], 
or within the framework of CDT. 

What is new in the CDT string field theory considered here is that we 
can use the corresponding stochastic field theory to solve the model, since 
we arrive at closed equations valid to all orders in the genus expansion. Let 
us translate equations (91) to /-space 



dG(l ,l;t) 
dt 

where the extended Hamiltonian 



-H(l )G(l ,l;t), (92) 



H ® = ~ l ^2 + Xl ~ 9sl 2 = H (l) - g s l 2 (93) 

now has an extra potential term coming from the inclusion of branching 
points compared to the Hamiltonian H (l) defined in (19). It is truely re- 
markable that all branching and joining is contained in this simple extra 
term. Formally H{1) is a well-defined Hermitian operator with respect to the 
measure (40) (we will discuss some subtleties in the next section). 

We can now write down the generalization of Wheeler-deWitt equation 
(40) for the disk amplitude 

H(l)W(l) = 0. (94) 
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Contrary to W\ UJ (/) appearing in (40), W{1) contains all branchings and all 
topology changes, and the solution is precisely (75)! This justifies the choice 
g s = 0/2 mentioned above. Recall that E = does not belong to the 
spectrum of H (l) since W (l) is not integrable at zero with respect to the 
measure (40). Exactly the same is true for the extended Hamiltonian H(l) 
and the corresponding Hartle-Hawking amplitude W(l) 
We have also as a generalization of (18) that 

G(l ,l;t) = (l\e~ tH ^\l ) (95) 

describes the nonperturbative propagation of a spatial loop of length Iq to a 
spatial loop of length / in proper (or stochastic) time t, now including the 
summation over all genera. 

9 The extended Hamiltonian 

In order to analyze the spectrum of H(l), it is convenient to put the differ- 
ential operator into standard form. After a change of variables 

I = l -z\ m = v^(z), (96) 
the eigenvalue equation becomes 

If = Ed>(z), H(z) = -\^- 2 + \\z 2 + A _ | z 4_ (9?) 

This shows that the potential is unbounded from below, but such that the 
eigenvalue spectrum is still discrete: whenever the potential is unbounded 
below with fall-off faster than — z 2 , the spectrum is discrete, reflecting the 
fact that the classical escape time to infinity is finite (see [52] for a detailed 
discussion relevant to the present situation). For small g s , there is a large 
barrier of height X 2 /(2g s ) separating the unbounded region for I > X/g s from 
the region < I < X/(2g s ) where the potential grows. This situation is 
perfectly suited to applying a standard WKB analysis. For energies less 
than \ 2 /(2g s ), the eigenfunctions (41) of H (l) will be good approximations 
to those of H(l). However, when / > X/g s the exponential fall-off of ipn\l) 
will be replaced by an oscillatory behaviour, with the wave function falling 
off only like l// 1 / 4 . The corresponding ip n (l) is still square-integrable since we 
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have to use the measure (40). For energies larger than A 2 / (2g s ), the solutions 
will be entirely oscillatory, but still square-integrable. 

Thus a somewhat drastic change has occurred in the quantum behaviour 
of the one- dimensional universe as a consequence of allowing topology changes. 
In the original, strictly causal quantum gravity model an eigenstate ipn\l) 
of the spatial universe had an average size of order 1 / \f\. However, allowing 
for branching and topology change, the average size of the universe is now 
infinite! 

As discussed in [52], Hamiltonians with unbounded potentials like (97) 
have a one-parameter family of selfadjoint extensions and we still have to 
choose one of those such that the spectrum of HQ) can be determined un- 
ambiguously. One way of doing this is to appeal again to stochastic quan- 
tization, following the strategy used by Greensite and Halpern [53], which 
was applied to the double-scaling limit of matrix models in [54, 55, 52]. The 
Hamiltonian (85) corresponding to the Fokker-Planck equation (91), namely, 

d V(^) d fdS(x) \ . . fx 3 \ , nn . 

is not Hermitian if we view x as an ordinary real variable and wave functions 
ip(x) as endowed with the standard scalar product on the real line. However, 
by a similarity transformation one can transform H(x) to a new operator 

H(x) = e- s(x)/29 °H(x) e s ^' 29s ] tjj(x) = e" 5(:c)/29 ^(x), (99) 

which is Hermitian on L 2 (R,dx). We have 

H{x) -- g '^ + {^{^) + 2^ i J' (1 ° 0) 
which after substitution of the explicit form of the action becomes 

H(x) = -g s ^- 2 + V(x), V(x) = ^-(A - x 2 ) 2 + x. (101) 
The fact that one can write 

R -^-i + ^rt (102 » 

implies that the spectrum of H(x) is positive, discrete and unambiguous. We 
conclude that the formalism of stochastic quantization has provided us with 
a nonperturbative definition of the CDT string field theory. 
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